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a . 

We study the geometrical optics generated by a refractive index of the form 
n(x, y) = 1/y (y > 0), where y is the coordinate of the vertical axis in an orthogonal 
■ reference frame in M 2 . We thus obtain what we call "hyperbolic geometrical op- 

tics" since the ray trajectories are geodesies in the Poincare-Lobachevsky half-plane 

in ■ 

H 2 . Then we prove that the constant phase surface are horocycles and obtain the 

00 

' horocyclic waves, which are closely related to the classical Poisson kernel and are 

the analogs of the Euclidean plane waves. By studying the transport equation in 
the Beltrami pseudosphere, we prove (i) the conservation of the flow in the entire 

x 

strip < y ^ 1 in H 2 , which is the limited region of physical interest where the ray 
trajectories lie; (ii) the nonuniform distribution of the density of trajectories: the 
rays are indeed focused toward the horizontal x axis, which is the boundary of H 2 . 
Finally the process of ray focusing and defocusing is analyzed in detail by means of 
the sine-Gordon equation. 
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I. INTRODUCTION 

It is well known that the geometrical optics approximation of the wave equation is related 
to the asymptotic form of the integral representation of the field (if such exists), which is 
an exact solution of the wave problem. Suppose, for instance, that the field in a uniform 
medium can be written as an expansion in plane waves; the evaluation of this integral by 
the stationary phase method yields an asymptotic series. Then the leading term of this 
asymptotic expansion, which is composed by an amplitude and a phase, can be extracted to 
yield the approximation. The ray trajectories are the lines orthogonal to the constant phase 
surface and are described by the eikonal equation; the amplitude satisfies the transport 
equation, whose physical meaning is related to the conservation of the flow. In the simplest 
case of uniform medium, whose refractive index n is a real constant, the rays are straight 
lines which are characterized by the following properties: 

(i) They are geodesies of the Euclidean space. 

(ii) Phase and amplitude are real-valued functions. 

(iii) They can be derived by the Fermat's principle. 

Constrained by these properties the methods of geometrical optics are rather limited and 
fail to explain several phenomena as, for instance, the diffraction by a compact and opaque 
obstacle, that is the existence of non-null field in the geometrical shadow which, for this 
reason, is usually referred to as the classically (or geometrically) forbidden region. 

In the decade 1950-1960 J. B. Keller-^ wrote several papers where he introduced the so- 
called Geometrical Theory of Diffraction (GTD). The latter can be regarded as an extension 
of geometrical optics, which accounts for diffraction by introducing the diffracted rays in 
addition to the usual rays of geometrical optics. After these seminal works there has been 
a steady flow of papers addressing various aspects of the theory. On the one hand papers 
oriented to pure and applied electromagnetic theory, like radiation and scattering of waves, 
antenna design, waveguide theory and so on- ; on the other hand, a highly theoretical and 
mathematically sophisticated theory of propagation of singularities and diffraction of waves 
on manifolds-. In spite of these efforts and a wide literature on these topics, not all the cases 
of interest have been studied. An example is what we could call the "hyperbolic geometrical 
optics" , that is the geometrical optics generated by the rays in the specific case of a refractive 
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index of the form n(x,y) — 1/y (y > 0), where y denotes a spatial coordinate, say vertical, 
in an appropriate orthogonal reference frame in IR 2 . As far as we know, this problem has 
never been treated, except for some very marginal remarks (see, for instance, Ref. la), in 
spite of its intrinsic geometrical interest and some possible applications to the physics of 
nonuniform optical fibers. It is precisely the main purpose of the present paper to fill this 



Let us return to Keller's program of widening, from a geometrical viewpoint, the arena 
of Euclidean geometrical optics. With this in mind we adopt, first of all, the Jacobi's form 
of the principle of least action (instead of Fermat's), which concerns with the path of the 
system point rather than with its time evolution^. More precisely, the Jacobi's principle 
(generally applied in mechanics) can be formulated as follows: If there are no forces acting 
on the body, then the system travels along the shortest path length in the configuration 
space. Here we assume a wide extension of Jacobi's principle, which can be formulated as 
follows: the geodesies associated with the Riemannian metric n(x, y) a/ dx 2 + dy 2 , i.e. the 
paths making the functional J n(x, y) \J dx 2 + dy 2 stationary, are nicknamed rays. In other 
words, in place of Fermat's principle which reads 5 Jp 1 dt = 0, where dt is the travel time 
measure, and Pq and Pi are prescribed starting and end points of the path, we write 



where y' = tana, a being the angle that the tangent to the curve y = y(x) forms with the 

X cLXIS. 

The simplest realization of this Jacobi's principle consists in identifying n 2 with the 
Riemann metric tensor g^, whenever this identification is admissible. This identification 
requires great caution, indeed; the form gijdx l dx^ must be symmetric and positive definite, 
and this poses a strict restriction. For instance, consider a refractive index (or, in mechanics, 
a potential) of the following form: n 2 = 1 — V/E, where E is the energy of the incoming 
particle and V is the height of the potential, with V > E as in the case of the tunnel effect. 
In this situation the geometric interpretation of the trajectory as a real-valued geodesic 
in a Riemannian manifold is no longer possible. The only chance remains to extend the 



gap. 




(1) 



or, equivalently, 
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admissible values of the phase to imaginary and/or complex values and, consequently, to 
speak of complex rays in the sense of Landau 7 . 

But let us return to the cases where this identification is admissible. As we already 
mentioned, it is obviously possible in the case of a uniform nonabsorbing medium: in this 
case we simply obtain a physical realization of Euclidean geometry. But it is also certainly 
admissible when the refractive index is of the form introduced above, i.e., n(x,y) = 1/y 
(y > 0), where y denotes the coordinate of the vertical axis in an orthogonal reference frame 
in M. 2 . In this case we are led to the Lobachevskian metric: ds 2 = (dx 2 + dy 2 )/y 2 . Then 
the rays are geodesies in the hyperbolic half-plane (Poincare half-plane): i.e., Euclidean 
half-circles with centers on the x axis (horizontal axis), or Euclidean straight lines normal 
to the x axis. Let us recall that the refractive index n is defined as: n = c/t> p h, where c is 
the light speed in vacuum, and v p h is the phase velocity of radiation of a specific frequency 
in a specific material. Therefore n ^ 1, and in the case n(x, y) = - only the strip < y ^ 1 
has physical interest; hence the actual rays will lie necessarily in this band. Accordingly, 
hereafter, the only optical paths considered will be the Euclidean half-circles with centers 
on the and radius R bounded by < R ^ 1. 

The subsequent step in developing an optical geometry consists in finding the constant 
phase surfaces and, accordingly, describing the analog of the Euclidean plane wave. This 
problem will be solved in Sec. [Til studying some geometrical properties of horocycles and 
introducing what we call horocyclic waves, which play in hyperbolic geometrical optics the 
same role as the plane waves do in the Euclidean one. At this point we have the main 
ingredients needed for writing the geometrical approximation of the wave function; what it 
is still missing is an analysis of the amplitude and of the related flux density. This latter 
problem can be analyzed at two different levels. First we prove that the flow of rays is 
conserved: once a pointlike source is fixed, no ray will be absorbed or created. This result 
will be proved in Sec. [TU A more subtle question is the following: Is the flow of the ray 
trajectories homogeneous or do the rays focus? This issue, besides its intrinsic geometrical 
interest, could in our opinion be of some interest in possible applications to the propagation 
in optical fibers with non-uniform refractive index-. This problem will be analyzed in detail 
in Sec. IIHI First we study the transport equation in the Beltrami pseudosphere, and prove 
that the flow of ray trajectories is not homogeneous, but there is a focusing of rays on 
the horizontal x axis. Glancing to possible applications to propagation in optical fibers 
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this result suggests a conjecture indicating a strong ray focusing along the fiber axis, when 
the refractive index profile in the fiber is of hyperbolic type, instead of paraboliclike, as is 
customary. Next, this problem will be reconsidered by studying the variation of the angle 
that the tangent to the meridian of the Beltrami pseudosphere makes with the rotation 
axis of this surface, which can be indeed represented as a surface of revolution generated 
by a curve in IR 3 . This leads to the sine-Gordon equation and provides a more precise 
description of the ray focusing and defocusing processes. This analysis is necessarily local, 
since the problem is worked out inside each horocycle; at the end of Sec. 1111} we show how 
to pass from a local description of the flow inside each horocycle to a global one. 

Finally, in the Appendix, the geometric and algebraic ingredients which occur in Sees. 
HT1 and II III will be given. This appendix is split in three parts: the first part is devoted 
to the various models of hyperbolic geometry and to the conformal maps which allows the 
transformation between them; in the second part we study the group SU(1, 1), which acts 
transitively on the non-Euclidean disk, and prove some relationships connecting the spherical 
functions to the horocyclic waves; the last part is devoted to the Beltrami pseudosphere. 



A. Variational minimization of the Jacobi's functional and the rays in hyperbolic 



Let us consider the upper half-plane model of the hyperbolic two-dimensional space H 2 : 
i.e., U = {z = x + iy : y > 0} equipped with the metric d derived from the differential 
ds = \dz\/lmz (see the Appendix). Then we apply the typical methods of variational 
calculus to the Jacobi functional 



(Pq and P\ denote two points of the ambient space where light propagates). First we prove 
the following proposition, which refers to the whole upper-half plane U. 



II. THE FLOW IN THE STRIP < y < 1 



geometrical optics 




or, equivalent ly, 
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Proposition 1. (i) Let J be the following functional 



Jx y 

and let F denote the integrand of |3|). The Euler-Lagrange equation for this functional reads 

yy" + y' 2 + 1 = 0. (4) 

(ii) The extremals of functional ([3J) are Euclidean half-circles with centers on the x axis, or 
Euclidean straight lines normal to the x axis lying in the half-plane y > 0. These are the 
geodesies in the hyperbolic geometry realized in the half-plane y > 0. 
(Hi) The Weierstrass condition for the functional ([3)] reads 

F yy = yi i + y2 )3 /2 > (v>°)> ( 5 ) 

and it is satisfied for any y' . 

(iv) There exists a field of extremals of functional (TJ|) ; and the transversality condition 
becomes an orthogonality condition of these extremals to the curve &(x,y) = const, (constant 
phase curve), which satisfies the following equation (eikonal equation): 



dxi dxj 



1, (6) 



where X\ = x, X2 = y, and gij is the metric tensor. 

Proof. The proof makes use of standard procedures and can be found, for instance, in Ref. 

B 

Remark. Let us recall once again that the domain of physical interest where the optical 
paths necessarily lie (in view of the fact that n ^ 1) is the strip < y ^ 1; therefore we 
shall consider only a subclass of the extremals of functional ([3]): i.e., the half-circles with 
centers on the x axis and radius bounded by < R ^ 1. 



B. Poisson kernel and horocyclic waves 

Let us now give a more precise formulation of the physical problem. Suppose that a 
pointlike source of light is pushed to — oo on the (x, y) plane. For the sake of simplicity, 
here we limit ourselves to the scalar representation of light, and phenomena associated with 
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FIG. 1: Horocyclic flow outside the forbidden region in the Poincare disk D. 



polarization will not be considered. From Proposition [T] it follows that light rays are half- 
circles with centers on the x axis. For several reasons which will appear clear in what follows, 
it is convenient to map conformally the half-plane y > into the unit disk \(\ < 1, which 
amounts to pass from the Poincare half-plane model U to the Poincare disk model D (see the 
Appendix and Fig. []]). The appropriate conformal mapping is given by: ( = i(z — i)/(z + i) 
(z = x + iy; ( = £ + irj). 

In the unit disk the light source will be located at ( = i. The band y > 1 will be mapped, 
in the (-plane, into the disk tangent to the boundary B of D in i with Euclidean radius |, 
and represents the forbidden region for the light rays. The circular arcs lying in the half- 
plane y > and normal to the x axis will be mapped, in the unit disk, into circular arcs 
perpendicular to the boundary |£| = 1, which are precisely the geodesies of the hyperbolic 
geometry in the unit disk model. 

From the transversality condition [see statement (iv) of Proposition [T] , it follows that the 
constant phase curve is the curve that intersects orthogonally the extremals of functional 
([3]): i.e., the geodesies. In the unit disk, parallel geodesies are geodesies corresponding to 
the same point b = e 1 * on the boundary B of D. Therefore, in the physical problem being 
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treated, the circles tangent to the unit circle at the point 6, which intersects orthogonally the 
pencil of parallel straight lines (i.e., arcs of circle orthogonal to B) are the constant phase 
curves, they family of horocycles, and are denoted by Hb- 
We can now state the following proposition. 



Proposition 2. (i) The Poisson kernel 

i-ICI 2 



P(C,b) 



(C=\C\e w ;b = e^), 



l + |C| 2 -2|C|cos(£- 
is constant on each horocycle Hb with normal b. 
(ii) The function 

i-ICI 2 



(?) 



[P{CM 



(8) 



1 + |C| 2 - 2|C| cos(6> - ■ 

is an eigenfunction of the Laplace- Beltrami operator on the hyperbolic disk D corresponding 
to the eigenvalue v{y — 1). 

(Hi) The hyperbolic waves (horocyclic waves) are represented by the following expression: 

i-ICI 2 



(9) 



Ll + |C| 2 -2|C|cos(£-0 
where ((,b) is the hyperbolic distance between the origin of D and the horocycle of normal b 
passing through ( G D. 

(iv) The conical functions J'_i +iA (coshr) (i.e., the first kind Legendre functions of index 
(— | + iA) (A 6 M)j can be represented by 



y_ i+iA (coshr) = / e <5-a)<C*> db (AgR,B = {( : |C| = 1}), 
2 Jb 



(10) 



and correspond to the fundamental series of the irreducible unitary representation of the 
group SU(1, 1), which acts transitively on the hyperbolic disk D. 
(v) The following equality holds: 



y_ i +iA (cosh r) = 7_ i _ iA (cosh r) (A G 



(11) 



Proof, (i) The level lines of the Poisson kernel P((, b) are the circles tangent to the unit 
circle at the point b = e 1 ^: i.e., the images of the horocycles Hb with normal b (see Ref. Isl). 
(ii) The Laplace-Beltrami operator Ad on the hyperbolic unit disk D is given byi^ 



a d = i [i - «^)r + £) 



(12) 
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If v G C is any complex number, a direct computation gives 10 

A d P v (C,6) = i/(i/-1)P"(C,6). (13) 

(iii) In the Euclidean case the function x — > e lk ^ ,U} \ where k G R, u> G S^ -1 ), x G M n , rep- 
resents a plane wave with normal us. It is indeed constant on each hyperplane perpendicular 
to uj, and furthermore is an eigenfunction of the Laplacian on M. n . The geometric analog 
of the plane wave in the case of the hyperbolic disk D is the function represented by the 
equality ([9]) (see Ref. Illl ). In fact, it is an eigenfunction of the Laplace-Beltrami operator 
on D, as proved by statement (ii) [see Eq. (1151) ]■ Further, putting 9 = in formulae ([7]) and 
flU), we have: 

in rT^Wj = ln i^H = rf(0 ' = (lc|e ^' = (c ' 6) ' (14) 

where d(0,() — hi[(l + |C|)/(1 — ICI)] is t ne hyperbolic distance between the origin and the 
point ( G D (see the Appendix). Therefore, ((, h) is the hyperbolic analog of (x, uj). In fact, 
in view of statement (i), ((, b) is the distance between the origin and the horocycle of normal b 
passing through £ G D, assuming that the origin falls outside the horocycle; ((, b) is positive 
if the origin is external to the horocycle, while it is negative (((, b) = ln[(l — |C|)/ (1 + ICI)]) 
if the origin is internal to the horocycle. 

(iv) If we put: £ = tanh(r/2) cos#, 77 = tanh(r/2) sin#, then \(\ = tanh(r/2). The Rieman- 
nian metric ds 2 = [4(d£ 2 + dr/ 2 )/(l — £ 2 — r/ 2 ) 2 ] becomes ds 2 = dr 2 + sinh 2 r d9 2 . By the use 
of this substitution in the expression of the Poisson kernel (JTj) or (jSJ), we have: 

1 — ICI 2 



l + |C| 2 -2|C|cos(#-0) 



[cosh r — sinh r cos( 



> G C), (15) 



and the integral sum of horocyclic waves [see statement (iii)] gives (see Ref. [ll| and Propo- 
sition [6] in the Appendix): 

I/ m<i " = T, C ( coahr + Ihroos H ** = T -^" Shr> <^ C >' < 16) 

where B is the boundary of the hyperbolic disk D, and ^^(coshr) are the first kind Leg- 
endre functions^. Finally, setting v = ~ — iA (A G R) we obtain the conical functions 
7_i +iX (cosh.r), which correspond to the fundamental series of the irreducible unitary rep- 
resentation of the group SU(1, 1): i.e., the group of the matrices of the forrn^ 



a c 
c a 
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\a\ 2 — \c\ 2 = 1; a, c G C, which acts as a group of isometries of the hyperbolic disk D by 
means of the map 

= (CeC) - (17) 

(v) Equality f TTTT) is proved in the Appendix (see Proposition [6]). □ 

Remark. It is well known that the classical Fourier transform refers to the decomposition 
of a function, belonging to an appropriate space, into exponentials of the form e lkx (k 
real), which can also be viewed as the irreducible unitary representation of the additive 
group of real numbers. Analogously, the exponentials e^ k,x ^ are characters of the group 
R 2 . But the hyperbolic disk is not a group. Therefore a straightforward generalization of 
the exponential for D is not possible. Nevertheless, in view of the fact that the function 
y_j,(coshr) corresponds to the fundamental series of the irreducible unitary representation 
of the group SU (1, 1) for v = \— iA, the exponential e^" 1 ^)^^ (A G R) represents the analog 
of the Euclidean exponential, and plays the same role in the hyperbolic Fourier analysis^. 

C. Conservation of the flow 

As already said in the Introduction, the ray trajectories are the lines orthogonal to the 
constant phase surface, and are described by the eikonal equation; moreover, {(, b) is the 
hyperbolic distance between the origin and the horocycle Hb of normal b passing through 
(. Therefore, in close analogy with the Euclidean optical geometry, and recalling that 
"y-i + i\{cosh.r) = !P_i_ iA (cosh r) (A G R) [see statement (v) of Proposition [2] , the expression 
of the analog of the Euclidean plane wave e lkx (k G R) can be written as follows: e^~ lA )^> fe ) 
(A G R). Thus the geometrical approximation of the wave function ip can be obtained by 
multiplying e^~ lA - ) ^'^ times a function which represents the amplitude. Then we can state 
the following proposition. 

Proposition 3. The geometrical approximation of the wave function ip reads: 

^(C,A,6) = A(\)e^- iX ^ (\ER,(£D,beB), (18) 
and the flow in the entire strip < y ^ 1 is conserved. 
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Proof. Let a be the conformal map 



Z = (7(C) 



• C + i 



(19) 



C-i' 



defined in the Appendix, that transfers the geometry of D into U . Since <j(0) = i and 
a(i) = oo, then the image by a of the horocycle Hi passing through ( = is the horizontal 
line Hoc = {x + iy : y — 1} in U (the horocycles in the Poincare half-plane will be 
hereafter denoted by H b ). The image by a of the horocycle H a -i^ tangent to Hi in D 



order to avoid proliferation of notations, we denote by the same letter b both the points on 
the boundary B of D and the corresponding points belonging to the boundary of H 2 , i.e. 
belonging to M). 

We already saw that the horocycle H b of normal b is perpendicular to each geodesic 
starting from b. To calculate the amplitude of the wave function, we must see how many 
geodesies perpendicular to H b intersect H b: with the additional condition that these geodesies 
belong to the band < y ^ 1. This corresponds to find the amount of normal vectors at 
H b , with unit norm, that are tangent vectors of geodesies in the band < y ^ 1. 

In general, if b is a point in RU {oo} and Till is the unit tangent bundle of U, then the 
horocycle flow hj :b : T-JJ — > T-JJ is the flow which slides the inward normal vectors to each 
H b to the right along H b at unit speed. To find the equation of the flow hj >b , first we consider 
the flow hj j00 of geodesies perpendicular to the horocycle ifoo of normal oo. Then we choose 
a transformation M b which maps the horocycle H^ into the horocycle H b . In particular, 
the map M b transfers the flow hj >oc into the flow hj >b . 
From the definition, 



where Vi denotes the unit vector vertically upwards based at i e U. This is because in the 
simplest case of horocycle flow hj jOQ , the geodesies perpendicular to H^ are vertical lines and 
the isometry sending one vertical line into another vertical line is the horizontal translation. 
Therefore, the horocycle flow along Hoc is simply the horizontal translation. 

Let us now consider the transformation M b such that M b (oo) = b. Then, the horocycle 



is the horocycle H b of radius | through b G K and tangent to the horizontal line ifoo (in 




(20) 
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flow hjb along H b is the image of hj tOC by M b , hence 

h jtb (v) = M b ( 1 J \v. (21) 




It is clear from the definition that the amount of geodesies in the flow hj )b does not depend 
on the radius of the horocycle. Given two different points b\ and 62 in the boundary of the 
hyperbolic plane, then the composition of M bl and M b ~ x sends the point 62 in &i- Moreover, 
M bl o M^ 1 sends the horocycle flow hj ib2 into the horocycle flow hj bl . This proves that the 
amplitude of the wave does not depend on b and £. 

Using Proposition [2l we obtain that there exists a function A(\) independent of ( and 
b such that Eq. (ITS]) is satisfied, and the conservation of the flow along the entire strip 
< y ^ 1 is proved. □ 



Remark. It is interesting to compare the propagation of light in vacuum with that within 
the strip < y ^ 1 belonging to H 2 . In vacuum each ray cuts orthogonally all the constant 
phase planes: i.e., each ray emerging from a plane cuts orthogonally all the other parallel 
planes. In H 2 propagation proceeds in a completely different form. Take two horocycles 
lying in the strip < y ^ 1, and tangent at the point z = (1 + i)/2: the first horocycle, 
denoted by Hq, has normal 60 = 0; the second one, denoted by Hi, has normal b\ = 1. Only 
one geodesic, denoted 7 t , lying in H , cuts orthogonally H\\ it emerges from b = and 
ends at 61 = 1. All the geodesies 7>, emerging from b = and lying in H above 74, cut 
orthogonally horocycles H b with b > 1; the geodesies 7<, emerging from bo = and lying in 
H below 7j, cut orthogonally horocycles H b with b < 1. However, the density of the flow 
of geodesies entering orthogonally each horocycle equals the density of the flow of geodesies 
exiting orthogonally the same horocycle. 

III. TRANSPORT EQUATION AND DISTRIBUTION OF THE DENSITY OF 

TRAJECTORIES 

A. Transport equation in the Beltrami pseudosphere 



Working out the problem in the space H 2 allows us to describe each trajectory as a 
geodesic in the Poincare plane (or disk), but this setting is not appropriate for describing 
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the evolution of a bunch of trajectories. Hereafter we will switch to a representation more 
suitable for an effective characterization of the amplitude factor in the geometrical approx- 
imation of the field. To this aim, let us first recall the following well-known negative result 
due to Hilbert: there is no regular smooth immersion X : H 2 — > M 3 . However, one can 
look for a local immersion X : U — > M 3 , where X is a continuous differentiable function, 
and U C i 2 is an open subset. We keep for U an open horocycle based at b. This local 
immersion can be realized by means of the Beltrami pseudosphere, denoted hereafter by Pj, 
(see the Appendix and Fig. [2]). In fact, let us consider in the hyperbolic disk D an infinite 
strip lying between two parallel straight lines emerging from the source point located on the 
absolute at C — — i- Then we take on these parallel geodesies a pair of points Aq and Bq, 
lying on a horocycle of normal 60 — e~ 17r / 2 = — i and cutting orthogonally these straight lines; 
Aq and Bq are spaced at distance of 2ir. One is then led to consider the domain (— i, A , Bq). 
The Beltrami surface cut along any of its generators can be isometrically mapped into the 
domain (— \,Aq,Bq) (see Ref. Il4l ). On a Lobachevskian plane there always exists reflection 
(i.e., a hyperbolic isometry) about an arbitrary straight line; in particular, reflecting the 
strip (— i, Aq, Bq) about the straight line (— i, Aq) we obtain a new strip isometric to the 
initial one and realized as a cut of the Beltrami surface in M. 3 . Reflecting then this new strip 
(— i, A\, Aq) (the segment AiA has length 27r) about the straight line (— i, A\) we obtain the 
strip (— i, A 2 , Ai) with the same properties. Exactly the same procedure can be repeated 
on the other side of (— i, Aq, Bq), leading to (— i, B2, B\). We thus obtain strips of the form 
(-i, A k , A k -i) and (-i, B k , B h _ x ) (1 ^ k < 00); all segments (A k , A k -i) and (B k , B k ^) have 
the same length 2ir. Working with the same procedure we can now construct the map of the 
open horocycle H bo , tangent at the boundary to the forbidden region (this latter represented 
by the horocycle if; of normal i and passing through the origin), into a Beltrami funnel, such 
that each strip of the type (—i,A k ,A k _i), (— i,A ,B ), and (—i,B k ,B k _i), (1 ^ k < 00) 
(referred, now, to the horocycle Hb Q ), is mapped isometrically into the Beltrami surface, the 
horocycle being wound infinitely many times into the Beltrami surface^ (see Fig. [2]). 
We can repeat the same procedure for each point b G B, since there is a rotation (i.e., a 
hyperbolic isometry) sending each b G B onto b . 



For an explicit equation of the immersion X, the reader is referred to Ref. 
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In general, the Laplace-Beltrami operator Am on a two-dimensional Riemannian mani- 



X 




FIG. 2: Mapping of a horocycle in the disk D into a Beltrami pseudosphere. 



fold M with metric tensor (g = |det(py)|, g u = ) is defined as follows: 



V 7 ^ 



• I 4 \j=l 



(22) 



In the specific case of the hyperbolic metric associated with the refractive index n(y) — 1/y 
(see the Appendix), the Laplace-Beltrami operator reads: 



A H = — 



1 / d 



d 2 



n 2 \dx 2 dy 2 
We then have the following proposition. 



d 2 



d 



32 



dx 2 dy 2 J 



(23) 



Proposition 4. (i) The Helmholtz equation reads 



A H ^ + kli> = 0, 



(24) 
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where = A 2 + \ (AG K). 

(wj The geometrical approximation of the wave function ip (for |A| — > oo^), written in terms 
of the Beltrami coordinates (see the Appendix), reads 

*fj±{\,u) = C{\)e u/2 e TiXu (AGl;u>0). (25) 

Proof, (i) Let us consider the horocyclic waves which generate the conical functions 
J , _i ±iA (coshr), corresponding to the irreducible unitary representation of the 577(1,1) 
group, which acts transitively on the hyperbolic disk D. This amounts to put in the expo- 
nent v E C of the Poisson kernel: v = |±iA (Ael). Accordingly, the horocyclic waves read 
e (^^ x )(CM [g ee statements (iv) and (v) of Proposition [2] . From statement (ii) of Proposition 
|2]and Eq. flU we get: 

AHe (i±iA)<C, fe > = _ ^2 + 1^ e (i ± iA){ Cl 6> = ^ e (|±iA){C, b > ) ( 26) 

where k^ = A 2 + \ (A G 1). Next, proceeding in close analogy with the Euclidean case, 
where the Euclidean plane wave plays the role of the horocyclic wave, we obtain Eq. (1241) . 
(ii) Let us now go back to the mapping of the horocycle into the Beltrami funnel (without 
a cut) in M 3 , illustrated above. Next, we apply the Laplace-Beltrami operator to the wave 
function ip, supposed to belong to C°°(S) (23 denoting the Beltrami pseudosphere); in (|22|) 
Xi (i = 1,2) stand for the Beltrami coordinates u,v. Recall that the first fundamental form 
in Beltrami coordinates reads [see part (C) of the Appendix]: 

/ = du 2 + e~ 2u dv 2 (u^O). (27) 

Accordingly, we have g n = 1, 5-22 = e~ 2u , g l2 = £21 = 0, g = \ det(^)| = e~ 2 ", = g^ 1 . 
Thus, we are led to the following equation: 

A 3 V + 4V> = 0, (28) 

where is the Laplace-Beltrami operator, referred to the Beltrami pseudosphere. In 
this equation, we pass from the coordinates (£, rf) of the hyperbolic disk D to the Beltrami 
coordinates (u, v) of the Beltrami pseudosphere. We illustrate with more details this passage. 
First we embed an open horocycle Hj, of normal b and tangent to the forbidden region 
(represented by the horocycle Hi passing through the origin of D and with normal i; see 
Fig. [1]) into a Beltrami pseudosphere. Notice that in the present analysis, as well as in 
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Proposition [31 and in strict analogy with the classical Euclidean procedure, we consider the 
distance from the origin of the hyperbolic disk D (rather than from the point source located 
at £ = i) to the horocycle H^ b (inside H b ) of normal b passing through a point (. Thus we 
have 

(C, b) = d(0, H b ) + d(H b , H Ctb ) := d b + d(H b , H c>b ). (29) 

When we embed H b into a Beltrami pseudosphere, the distance d(H b , H^ b ) between horo- 
cycles corresponds to the distance between different parallels u = const, inside the pseu- 
dosphere. Since H b is fixed, then d b is fixed too. Then, following the standard method of 
stationary phase, we look now for a solution of equation (128]) . of the following form: 

V>(A,x) = [ A(x,£)e^- im ' b) d£ = e^-' lX)db [ A(x, i) e (^)d(H b ,H Cb ) d£ 



= C(A) J A(x,£)e^- iX ^ x ' e) d£, (30) 

where x = (xi, x 2 ), x\ = u, x 2 = v ; £ is the pathlength inside the pseudosphere to the point 
of coordinate x, $(x, £) denotes the phase (recall the statements of Proposition [2]). 

The r.h.s. of equality (1301) is an integral of oscillating type. The principal contribution 
to x), as |A| — > +oo, corresponds to the stationary point of $, in the neighborhood of 
which the exponential ceases to oscillate rapidly. These stationary points can be obtained 
from the equation d§/d£ = (provided that d 2 $/d£ 2 ^ 0). If the condition d§/d£ = 
is satisfied by a unique value £q of £, corresponding to the unique ray trajectory (geodesic) 
passing across the point of coordinates (u,v), we say that $ has a critical nondegenerate 
point at £ = £q. Moreover, recalling that the manifolds with nonpositive curvature do not 
have conjugate points, we can state that all the critical points of $ are nondegenerate. Then, 
by applying the Morse lemma on the representation of the functions all of whose critical 
points are nondegenerate, we obtain the following asymptotic evaluation of integral ( |30l) : 



V>(A, x) = C(A)e**<*A> e~™M £ (31) 
The leading term of expansion (13T|) reads 



^(A, x) = C(A)A (x)e^ (xA) e- iA$(xA) , (32) 

where 



A (x) = A(x,4) 



d 2 $ 



d£ 2 



-l/2> 



cxp 



7T fd 2 § 



i — sgn 



4 & \d£ 2 



(33) 
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For simplicity, in the following we shall write the leading term of expansion (!3T|) as: 
C(A)A(x)e(2- iA )' I> W, dropping the zero subscripts. Substituting this expression into Eq. 
( 1281) . collecting powers of (iA) and, finally, equating to zero their coefficients, two equations 
are obtained: the eikonal (or Hamilton- Jacobi) equation 



and the transport equation 



J_ V — 



dxj 



i=i J . 



(34) 



0. 



(35) 



Let us note that in the present problem the wave functions are radial, in view of the fact 
that we are considering a family of horocycles having all the same normal b. Therefore, 
■0(A, x) [where x = (u, v), u,v being the Beltrami coordinates, see the Appendix] does not 
depend on v. Then, Eq. f l34|) becomes: 

(36) 



d S i =1. 



du 

which gives the following expression of the phase: <3>( ± - ) = ±u + c (c = const.). Proceeding 
analogously with Eq. (135]) . we have: 



d 

du 



0. 



(37) 



Substituting in the leading term 
and ([37]), the rhs of Q25D follows. 



the expressions of $ and A, which derive from fl3T)|) 

□ 



B. Sine Gordon equation and the flow of trajectories 

The analysis of propagation in the Beltrami pseudosphere has allowed us to study the 
distribution of the density of trajectories and, accordingly, the ray focusing along the hor- 
izontal axis. Another parameter, related once again to the Beltrami pseudosphere, and 
whose characterization is relevant in our description of the flow of trajectories, is the angle 
ip that the tangent to the meridian (in the Beltrami pseudosphere) makes with the z axis 
[see section (C) of the Appendix]. 
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In the Appendix, the Beltrami pseudosphere is described in terms of the Beltrami coor- 
dinates u and v . Here we choose another parameterization (p, q) by setting 



Integrating (I38al) . we obtain: 



and 



dp 
dq 

¥ = 



= — esc ip dip, 
= —dv. 

2 tan- 1 (e~ p ). 



dip 

— = — sin ip 
dp 



coshp 



(38a) 
(38b) 

(39) 
(40) 



Next, substituting in formula (1A.31I) of the Appendix (with p = —1) p> and v in terms of the 
parameters p and q, the position vector r of the pseudosphere can be rewritten as follows: 



/ 



r(p, q) 



coshp 
1 



- — cos q ^ 



, sin q 

cosnp ^ 

, p — tanhp ; 



(41) 



where the downward vertex of the pseudosphere corresponds to p — > — oo, while the rim 
corresponds to p — 0. The parameters (p,q) can be related to the arc lengths (cx,j3) along 
asymptotic lines as follows: 

p = a + (3, 
q = a — (3. 

Setting lo = 2tp, the first fundamental form reads 



(42) 



I = da 2 + 2 coscuda d/3 + d/3 2 



(43) 



which can be easily derived from the expression (IA.32I) of the Appendix through the formulae 
(I38al) and (|42|) (recalling that p 2 = 1). Moreover, from Gauss' and Weingarten's equations, 
it follows that u, which is the angle between the asymptotic lines, satisfies the classical 
sine-Gordon equation, which reads 

d 2 u 



dad (3 smuJ - 

This latter equation rewritten in terms of the coordinates p, q becomes 

d 2 (p d 2 p> 

We can finally state the following proposition: 



(44) 



(45) 
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Proposition 5. (i) The angle ip, that the tangent to the meridian of the pseudosphere makes 
with the z axis (i.e., the rotation axis), is represented by the following formula: 

if = 2 tan" 1 (e~ p ), (46) 

which is the so-called "one-soliton" solution of the sine-Gordon equation d45j) . 

(ii) The angle (p varies from tp = | to ip = it in the process of focusing, and from tp — it to 

ip = | in the process of defocusing. 

Proof, (i) The proof of this statement follows easily by direct calculation, 
(ii) The downward vertical z axis of the pseudosphere we are considering is negatively 
oriented; then varying p from to — oo [in equation fj^BT) ]. tp varies from (p = | to if = tt 
(focusing). Next varying p from — oo to 0, <p varies from tp = n to tp = | (defocusing). □ 

So far in this section we have considered a local description of the flow in order to show 
that the density of the flow is not homogeneously distributed. We now want to recover a 
global description in the entire strip < y ^ 1. But we have already remarked that the 
map of a horocycle into a pseudosphere excludes the boundary of the horocycle. Therefore 
this global description cannot be reached by considering only horocycles lying in the strip 
< y ^ 1 and tangent to the line y — 1. Recall, however, that while the physical geodesies, 
which lie within the strip < y ^ 1 in HI 2 , cannot enter the forbidden region y > 1, this 
is not the case for horocycles, provided we limit ourselves to consider in these horocycles 
(i.e., those entering the forbidden region) those segments of geodesies which lie in the strip 
< y ^ 1. In view of these considerations it is sufficient to consider horocycles H' b D (H' b 
entering the forbidden region), and, accordingly, the maps X' b embedding H' b into a Beltrami 
pseudosphere P b '. We can thus obtain a sequence of Beltrami pseudospheres P b which, by 
varying 6, allows us to connect in M 3 the geodesies lying in the strip < y ^ 1. We thus 
pass from a local to a global description. 

IV. CONCLUSIONS AND DISCUSSION 

In this paper we have presented the geometrical optics generated by a refractive index of 
hyperbolic type. The ray trajectories are geodesies in the Poincare-Lobachevsky half-plane, 
and the horocyclic waves, which are related to the Poisson kernel, represent the analogs of 
the Euclidean plane waves. We thus obtain two main results: 
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(a) The flow in the entire strip < y ^ 1 is conserved (see Proposition [3]); 

(b) inside each horocycle (by embedding horocycles in Beltrami pseudospheres) the ray 
focusing on each point b of the horizontal x axis: i.e., toward the boundary of H 2 , is 
shown. 

The connection between hyperbolic geometry and optics of spatially nonuniform media 
can be made even tighter. In fact, it can be shown that the transfer matrix associated with 
lossless layered optical media is an element of the group SU(1, 1), no matter how complicated 
the stepwise profile of the refractive index might be^^^-. Therefore, the action of any 
lossless optical multilayer can be regarded as a Mobius transformation on the unit disk—, 
and therefore the natural geometric environment for these physical systems is the hyperbolic 
one. From this point of view, the geometrical optics description of light propagation in the 
"hyperbolic glass" discussed so far can be regarded as the study (in a spatially continuum 
setting) of the particular case of special interest, in which only motions along hyperbolic 
geodesies are allowed. 



Appendix 

(A) Let us consider the upper half-plane model of the hyperbolic two-dimensional space 
U = {z = x + iy : y > 0}. Then the boundary dll of U is the real axis and infinity. On U we 
can define a metric d derived from the differential ds = \dz\/lmz where dz is the standard 
Euclidean metric. The geodesies of U are vertical half-lines and Euclidean semicircles with 
center on the real axis. The group of the orientation preserving isometries of (U, d) is the 
Mobius group PS'Z^M), that is the group of the 2x2 matrices of real coefficients with 
determinant 1. The action on U UdU is defined as j(z) = (az + b)/(cz + d), where ( ^ ^ 
G PSL 2 (R). 

Another model of the hyperbolic two-dimensional space is the Poincare disk D = = 
£ + irj : £ 2 + rf < 1}. The map 

C = i^, (A.l) 
z + 1 

{with inverse z = — i[(C + i)/(C — i)]} transfers the geometry of U into the geometry of D. 
In particular, the metric on D is given by the differential 

d ^ = T^W dEC (A ' 2) 
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where d^ is the standard Euclidean metric. In this model, the geodesies are circular arcs 
perpendicular to the boundary \(\ = 1. Moreover, the group of the orientation preserving 
isometries of the Poincare disk D is the group 577(1, 1) of the maps of the following form: 

° C + S - 'a| 2 -|c| 2 = ll. (A.3) 



c( + a 

(B) The group G = SU(1, 1) admits two subgroups relevant for our analysis: 
(1) The subgroup K of rotations 

/ oif n \ 

(0 < 9 < 4tt). 




(2) The subgroup A of matrices 

cosh £ sinh ^ 

a,- = I I (/' e 

sinh | cosh ^ 

Let v4 + denote the set a r with r ^ 0. Then the following decomposition holds. 
Cartan decomposition 10 : Any element g G G can be decomposed as follows: 

g = k e a r k^ (r ^ 0, ^ 6 < 4tt, ^ < 2tt), (A.4) 

i.e., G = The decomposition is unique if g G" K. 

We now introduce the so-called spherical functions $ v {g) on G/K [g e G = 577(1, 1), f G 
C], which are defined as follows^. 

Definition 1. The spherical functions on G/K [G = SU(1, 1), K = 5*0(2)] are defined as 
follows: 



B 



dig- 1 ■ b) 



db [g G SU(1, 1), v G C], (A.5) 



db 

where B is the boundary of the hyperbolic disk D (i.e., B = {( : \(\ = 1} ). 



We can then prove the following proposition. 

Proposition 6. The functions $v(g) satisfy the following properties: 
(i) 

**(g) = ^(coshr), [g G SU(1, 1), v G C], (A.6) 
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where 7- v {-) are the first kind Legendre functions, 
(ii) 

A D <f>„(g) = u(u - l)$ v (g), [g e SU(l, l), v e 

where Ad is the hyperbolic Laplace- Beltrami operator. 
(Hi) 

y_i +iA (coshr) = y_i_ iA (coshr) (A G 1). 



(A.7) 



(A.8) 



Proof. Let us consider the following integral j- J Q 27r /(e 1 ^) d(p [/ G L X (B)], and evaluate how 
the Lebesgue measure ^d0 changes when an element g G SU(1, 1) acts on B (boundary of 
D). Recalling that the action of g is g • £ = (a( + c)/(c( + a) (£ G -D) (see Proposition [2]), 
we have g ■ e 1 ^ = e lx = (ae 1 * + c) / (ce 1 * + a) . We thus have^: 



dx 



ce lx — a\ 



—e 



C : 



-e 



(A.9) 



since a 



1. Let us now note that (7 • = c/a; on the other hand, in view of the 



Cartan decomposition, we have: 



g ■ = kg a r ■ = e ld tanh 



Therefore (1 - \c\ 2 /\a\ 2 ) = 1 - |C| 2 , and 

c ■ " 2 1 
-e 1 * - 1 
a 



IC|e i( 



1 



(A.10) 



1 
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2- 



ix 1 



d0 



|l_| C | e i(x-e)| 2 1 + |C| 2 -2|C|cos( X -6>)' 
We thus have 

and, in view of formulae ([Oil . (IQOjl . (jXTTjl . (tQ2j> . we get 

I-ICI 2 d(^-6) 



dx 



P(<?-0, 6) 



l + |C| 2 -2|C|cos( X -^) 
Recalling Definition [T] we can write 

1-ICI 2 



db 



(A.ll) 



(A.12) 



(A.13) 



L 


dig- 1 ■ b) 


v 1 p2ir 

db = — 




db 


2W 



l + |C| 2 -2|C|cos, 



d(f>. 



(A.14) 



Finally, writing tanh(r/2) in place of |£| [see formula flA.lOp ]. we obtain 

«2- 



1 



(cosh r + sinh r cos 1 



d(p = T_^(coshr), 



(A.15) 
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where the last equality follows from the integral representation of the first kind Legendre 
functions 12 . Formula (I A. 151) proves statement (i). 

From formula (1 A. 151) it follows that the spherical functions <& v {g) are bi-fc-invariant; 
indeed, using the Cartan decomposition, we have 

$ v (g) = ^ v {k 6 a r k^) = <$> u (a r ) = IP-^coshr). (A.16) 

Furthermore, we can prove that [see formulae ( 1A.13I) and ( 1A.14I) ]: 

$v(E)= [ [P(E-0,b)] v db = 1 = 3>_„(1). (A.17) 

J B 

Let us now return to statement (ii) of Proposition [2} and recall that -P"(£,6) is an eigen- 
function of the hyperbolic Laplace-Beltrami operator with eigenvalue v{y — 1). We then 
consider an integral of the following form: 

/ [P(g-0,b)] v db = $ v (g). (A.18) 

J B 

It can be proved that this integral superposition is still an eigenfunction of the hyperbolic 
Laplace-Beltrami operator A^, having v{y — 1) as eigenvalue^. In particular, it follows 
that § u (g) and [g G 577(1, 1)] are both eigenfunctions of A D with the same eigen- 

value, and therefore they coincide. We can thus state that T_i +iA (cosh r) = y_i_ iA (coshr); 
statements (ii) and (iii) are thus proved. 

Finally, from the representation of the Poisson kernel in terms of horocyclic waves, we 
have 

y_i , iA (coshr) = — / — d0 

2 Z7r J \ cosh r + smh r cos <p J 

2tt/ I/-I2 \ 1/2-iA 



— f ( , 1 ^) 2 , 1 d0= [ e^- iX) ^ b) db (A G R, C e D).(A.19) 

2vr7o Vl + KI 2 -2|C|cos0; * J B 1 



□ 



(C) In Sec. Hilt we have seen that it is possible to have a local isometric immersion from an 
horocycle to M 3 . The image of this immersion is a Beltrami pseudosphere. Now, we describe 
in more details the equation of the pseudosphere and its coordinates. 

First, we want to recover the pseudosphere as a surface of revolution of a curve in M 3 . 
Let us recall that, in general, the position vector r of the surface of revolution generated by 
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the rotation of a plane curve z = $(r) about the z axis is given by: 

' r cos v ^ 

r sin f , (A.20) 

where v varies between and 2ir. Here the circles r = const, are the parallels and the curves 
v = const, are the meridians. The first fundamental form associated with the surface (1A.20j) 
is given by 

/= {1 + [$'(r)] 2 }dr 2 + rW. (A.21) 
We now rewrite the form (IA.21I) as follows: 



I = du 2 + r 2 dv 2 , 



(A.22) 



where 



du = v 7 ! + [$'(r)] 2 dr; r = r{u). (A.23) 
From the general Gauss' theory of surfaces, we have that the total curvature is given by 20 



K 



1 d 2 r 



r du 2 



(A.24) 



whence the general pseudospherical surface of revolution with K = — 1/p 2 adopts the form 2 ^ 

(A.25) 



/ \ . u . _ u 

r(u) = c\ cosh — h C2 sinn — . 

P P 



In the case c\ = C2 = c, which corresponds to a parabolic pseudospherical surface of revolu- 
tion, the meridians are given by 

r( u ) = ce u/p , (A.26) 



while 



1-1-1 e 2u /p du. 



z = $(r) = 

Then, the first fundamental form, with c = 1, has the following expression: 

I = du 2 + e 2u/p dv 2 . 

The coordinates u, v are called Beltrami coordinates. 
The substitution 



(A.27) 



(A.28) 



sin ip = —e 
P 



u/p 



(A.29) 
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in f[A27l) yields 



z = p (cos <p + In 



tan 



(A.30) 



From formulae f[Q0|) . f[P6l) . fTAT29]) and fTA~30]) we obtain 



/ 



p sin <p COS V 
p sin <p sin v 
y p (cos<p + In |tan f |) y 

and the first fundamental form (in terms of <p and v ) is 



(A.31) 



I = p 2 cot z (p c\ip A + p 2 sin 2 <p dv 



(A.32) 



Equation (1A.31[) is the parametric form of the parabolic pseudosphere, seen as surface of 
revolution about the z axis of the curve called tractrix, which satisfies the following property: 
the length of the tangent from the point where it touches the curve to the point where it 
intersects the z axis is constant and equal to |p|; <p is the angle that the tangent to the 
meridian makes with the z axis. The angle <p varies between and tt so that, keeping p = 1, 
the related parabolic pseudosphere has vertices at z — +oo (corresponding to <p = n) and at 
z = — oo (corresponding to <p = 0) and rim at z = (corresponding to tp — |). The curve is 
continuous and regular except at the point z = 0, which is a cusp point. Choosing p = — 1 
and varying ip from to 7r, we shall have the upward vertical z axis positively oriented 
(<p varying from to |), and the downward vertical z axis negatively oriented (<p varying 
from | to 7r). In accordance with Hilbert's theorem, for which it is impossible to embed 
the entire hyperbolic disk onto R 3 , and since we want the immersion from the horocycle to 
the pseudosphere to be regular, then the image must be contained either in the downward 
component or in the upward component of the pseudosphere; thus it does not contain the 
cuspidal rim. Accordingly, taking p = — 1 once and for all, the first form, written in Beltrami 
coordinates, reads 

I = du 2 + (e~ u dv) 2 (u > 0), (A.33) 



(A.34) 



and the tractrix is: 

x = — sin ip, 

z = — (cos<p + In | tan ||) . 
Then, in Sec. IIII Al we use the form ( 1A.33j) with u ^ 0, varying ip from | to and, 
accordingly, z ^ (see Fig. [2]). In Sec. IIII where the coordinate u does not enter the 



26 



game and we have chosen the coordinates p and q, it is convenient to vary ip from | to n 
and, in accordance, taking the downward vertical axis negatively oriented, i.e., z ^ 0. 
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